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Neutron and x-ray powder and single crystal synchrotron diffraction of CsyFe2-xSe2 show the 
presence of superstructure reflections with propagation vector k=[|, |, 1] with respect to the average 
crystal structure M/mmm (a = 4, c = ISA). The propagation vector star corresponds to the 5 times 
bigger unit cell given by transformation A=2a-|-b, B= -a+2b, C= c. A solution for the atomic 
structure is found in the space groups P42/n and 14/ m with an ordered pattern of iron vacancies 
corresponding to the iron deficiency x — 0.29 and Cs stoichiometry y = 0.83. The superstructure 
sateflites are more pronounced in the neutron diffraction patterns suggesting that they can have some 
magnetic contribution. We have sorted out possible symmetry adapted magnetic configurations and 
found that the presence of AFM ordering with the ordered magnetic moment of Fe with ~ 2(Ub does 
not contradict to the experimental data. However, the solutions space is highly degenerate and we 
cannot choose a specific solution. Instead we propose possible magnetic configurations with the Fe 
magnetic moments in (a6)-plane or along c-axis. The superstructure is destroyed above ~ 500 K 
by a first-order-like transition. 

PACS numbers: 75.50.Ee, 75.25.-j, 61.05.C-, 74.90. +n 



I. INTRODUCTION 

The recent discovery of the Fe-based superconductors 
has triggered a remarkable renewed interest for possible 
new routes leading to high-temperature superconductiv- 
ity. As observed in the cuprates, the iron-based supercon- 
ductors exhibit interplay between magnetism and super- 
conductivity suggesting the possible occurrence of uncon- 
ventional superconducting states. Other common prop- 
erties are the layered structure and the low carrier den- 
sity. Among the iron-based superconductors FeSe has 
the simplest structure with layers in which Fe cations 
are tetrahedrally coordinated by S^. Recently super- 
conductivity at about 30K was found in XyFe2-xSe2 for 
X=K, Cs, RtPS. Muon-spin rotation/relaxation (/iSR) 
experiments evidence that the superconducting state ob- 
served in CsyFe2_xSe2 below 28.5(2) K is microscop- 
ically coexisting with a magnetic phase with a transi- 
tion temperature at T„i = 478.5(3) KP. The magnetic 
phase appears characterized by rather large static iron- 
moments as the fiSK signal is wiped out below Tm- Very 
recently the AFM order was reported- ' in superconduct- 
ing K0.8Fe1.6Se2 with Tn = 560 K with the iron magnetic 
moment 3.31 /is. 

The average crystal structure of XyFe2-xSe2 is the 
same as in the layered (122-type) iron pnictides with 
the space group lA/mmrn^. Different types of iron va- 



cancy ordering in TlyFe2-xSe2 were observed long time 

agcP^, including the one with 5 times bigger unit cell. 

Due to renewed interest to the superconducting chalco- 

genides many new experimental studies on the vacancy 

orderinginXvFe2-xSe2 (X=K,T1) have appeared very 
recentljiSlMiijiimil. 

In the present paper we report on the observation 
of superstructure in superconducting (Tc = 28.5 K) 
CsyFe2-xSe2 below Tg ~ 500 K and analyze the diffrac- 
tion data assuming iron vacancy ordering and possible 
antiferromagnetic ordering of Fe at room temperature. 
The single crystals used in the present study are the same 
as in the Refs.-^ 



II. SAMPLES. EXPERIMENTAL 

Single crystals of cesium intercalated iron selenides of 
nominal compositions Cso.8(FeSeo.98)2 were grown from 
the melt using the Bridgman method as described in 
Ref.'^. Powder x-ray diffraction was performed using a D8 
Advance Bruker AXS diffractometer with CuKq radia- 
tion. For these measurements a fraction of the crystal was 
cleaved, powdered, and loaded into the low background 
airtight specimen holder in a He-glove box to protect the 
powder from oxidation. Differential scanning calorimetry 
(DSC) experiments were performed with a Netzsch DSC 
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204F1 system. Measurements were performed on heating 
and cooling with a rate of 10 K/min using 20 mg samples 
encapsulated in standard Al crucibles. An argon stream 
was used during the whole experiment as protecting gas. 
Neutron powder diffraction experiments were carried out 
at the SINQ spallation source of Paul Scherrer Institute 
(Switzerland) using the high-resolution diffractometer for 
thermal neutrons HRPT"-^ (A = 1.866, 1.494 A, high in- 
tensity mode Ad/d > 1.8 • lO"'^). Refinement of crystal 
and magnetic structures of powder neutron diffraction 
data were carried out with FULLPROfE^ program, with 
the use of its internal tables for scattering lengths and 
magnetic form factors. Single crystal diffraction data 
were collected at the SNBL beamhne BMIA at the ESRF 
synchrotron in Grenoble (France) with a MAR345 image- 
plate area detector using A = 0.6977(1) A. Intensities 
were indexed and integrated with CrysAli empirical 
absorption correction was made with SADAB^, struc- 
ture refinement with SHELXL97iS-. 



III. RESULTS AND DISCUSSION 
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FIG. 1: Fragment of the neutron diffraction pattern of 
CsyFe2-xSe2 . Contribution of the superstructure peaks re- 
fined in the powder profile matching mode for k=[|,|,l] 
with respect to the average crystal structure {M/mmm with 
a = 3.96, c = 15.29 A and the structure parameters fixed by 
the values from Table [l| is shown by red solid curve. 

The average crystal structure can be refined in the 
standard structure modeP. The iron site occupancy is 
refined to smaller than unity for both x-ray and neu- 
tron diffraction data indicating the presence of the va- 
cancies on the iron sites. The structure parameters re- 
fined in this model are presented in Table |l] for both 
laboratory x-ray and NPD data. The neutron diffrac- 
tion pattern has a set of extra diffraction peaks that 



TABLE I: Crystal structure parameters refined in the av- 
erage structure model lA/mmm (no. 139), Fe in (0, |, |) 
(4d), Se in (0,0,^) (4e) and Cs in (0,0,0) (2a) positions. 
In the space group I4/m (no. 87) which unit cell is gener- 
ated by the transformation given in the text the atoms are 
split in the following way: Cs in (0, 0, 0) (2a) and (0.4, 0.8, 0) 
(8h), Se in (0.4,0.8,-2se) (16i) and {\,\,-zse + \) (4e), 
Fe in (0.3,0.6,0.25) (16i) and (|,0, i) (4d), where 0se is z- 
component in M/mmm. The stoichiometries o-Cs and o-Se 
are calculated to be in units of the formula CsyFe2-xSe2 . 
The data for laboratory x-ray and neutron powder diffraction 
NPD are given at room temperature, for the synchrotron sin- 
gle crystal (s.c.) experiment at T=536 K in true M/mmm 
symmetry above the order-disorder transition. 





x-ray 


NPD 


NPD,/4/m 


s.c. x-ray 




300 K 


300 K 


300 K 


536 K 
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3.9608(2) 3.9614(2) 


8.8582(3) 


4.0177(5) 
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15.285(1) 


15.2873(9) 15.2873(9) 


15.333(4) 


z-Se 


0.3456(4) 0.3436(3) 


-0.343(3) 


0.3443(3) 


o-Cs 


0.636(13) 0.622(24) 


0.73(1) 


0.754(9) 


o-Fel 


1.48(3) 


1.49(3) 


1.50(2) 


1.66(8) 


o-Fe2 






0.024(8) 




B-Fe 
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1.3(2) 


B-Se 
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3.2(1) 
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5.9(5) 


3.2(2) 
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5.95 


7.47 


7.12 
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x' 


3.6 


3.3 


2.86 


1.4 



can be indexed with the propagation k- vector k=[|, |, 1], 
as shown in Fig. [T] However the x-ray powder diffrac- 
tion pattern contains only one clearly visible satellite 
(|,— i,0) at q = 1 A^^, that allows one to suggest 
that the satellites seen by neutron might have mag- 
netic contribution. The propagation vector star k = 
{[f,i,±l],[-i,-i,±l],[i,-i,±l],[-i,|,±l]} corre- 
sponds to the new unit cell given by the transformation 
A=2a+b, B= -a+2b, C= c (the illustration of the lat- 
tice cell transformation is shown in Fig. [2]). A good refine- 
ment of NPD pattern explaining the satellite peaks can 
be done with the supercell indicated above in the space 
groups P4:2/n and 14,/ m. Using the fixed new atomic 
positions generated from the average crystal structure 
(lA/mmm) by applying the above basis transformation 
and releasing only the site occupancies and z-Se, simi- 
lar as for the average structure, one immediately gets a 
reasonably good description of the superstructure peaks. 
Table |l] shows the atomic positions and details of the re- 
finements. The Fe site splits in two sites in lA/m (no. 
87) and in three sites in P4:2/n [no. 86, note to work in 
the second setting with origin at —1 post-matrix transla- 
tion (i, i, i) should be applied]. Both groups give sim- 
ilar quality of the refinements of the single crystal data 
as we explain below, so we present the results only for 
lA/m space group, which is more symmetric with respect 
to the iron sites. The atomic positions for lA/m space 
group generated by the above transformation from aver- 
age space group lA/mmm are listed in Table |l] All sym- 
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FIG. 2: The lattice unit cell transformation to a new bigger 
tetragonal supercell shown by capital letters. The supercell 
shown by solid red lines corresponds to the propagation vec- 
tor star generated by ki = [|, |, 1]. The dashed cell shows 
the twin domain that corresponds to the star generated by 
k2 = [|, |, 1] The k- vector stars are shown in Fig. [s] Iron 
vacancy ordering pattern in ab-plane is shown by blue circles. 
The brown and green circles show fully occupied Fe and Se 
positions projected to the afe-plane. 




FIG. 3: A slice of the reciprocal space showing [hkO] plane. 
The indexing is given in the average cell (M/mmm). The 
satellite reflections are indicated by red arrows for ki = 
[|, |, 1] and k2 — [|, |, 1] by solid and dashed lines, respec- 
tively. These two k-vectors correspond to two twin domains 
shown in Fig. |2] 



metric" positional parameters which were generated from 
the special coordinates of lA/mmm were fixed in the re- 
finement. The fully occupied Fel (16i) site contributes 
1.6 of iron stoichiometry in CsyFe2-xSe2 . If both site 
occupancies are refined, Fel gets almost maximal value, 
whereas Fe2 (4d) site occupancy is close to zero, as shown 
in the Table |ll The iron vacancy pattern looks like as 
shown in Fig. [2] 

Since both vacancy superstructure and possible mag- 



TABLE II: Crystal structure parameters refined in the space 
group 74/m (no. 87) from the synchrotron single crystal ex- 
periment at room temperature. The Wyckoff site symmetry 
positions are the same as described in the caption of Table [l] 
The refined stoichiometry is Cso.g3(i)Fei.7i(i)Se2, The Fel site 
occupancy was fixed to 1. Anisotropic atomic displacement 
parameters Uij are in A multiplied by 10'^. U^q is defined as 
one third of the trace of the orthogonalized Uij tensor. To- 
tally 40 parameters were refined using 6286 refiections 539 
of which are independent. Final R factors are Rl = 0.0860, 
wR2 = 0.1955 [/ > 2o-(7)] and Rl = 0.0940, wR2 = 0.2105 
(all data). 

X y z Uiso occ 

Csl 0.0000 0.0000 0.0000 78(2) 0.911(14) 

Cs2 0.4041(2) 0.8057(2) 0.000 75(2) 0.81(1) 
Sel 0.3924(2) 0.7987(2) 0.6551(2) 55(2) 1 
Se2 0.5000 0.5000 0.1488(2) 54(2) 1 

Fel 0.3014(1) 0.5938(1) 0.25165(6) 59(2) 1 
Fe2 0.5000 0.0000 0.2500 56(9) 0.27(2) 

C^ii U22 Uz'i U23 U13 U12 
Csl 92(2) 92(2) 49(4) 
Cs2 84(2) 93(2) 49(3) 3.9(7) 
Sel 61(2) 59(2) 45(2) 2.1(5) 1.1(5) 1.6(5) 
Se2 58(2) 58(2) 44(3) 
Fel 63(3) 62(2) 51(4) 0.8(7) -1.1(6) -0.1(7) 
Fe2 62(10) 62(10) 45(16) 



TABLE III: Irreps Ti,i — 1..8 of the space group M/m (no. 
87) for k-vector k = 0. Numeration of the irreps is in accor- 
dance with Kovalev book (fc — kn, table ri21)^^. T4, = T3T2, 
Te = T5T2, Ts = T7T2. The basis functions ip for the spin of iron 
atom in general (16i)-position for T4 (two in (afe)-plane) and 
for T7 (one along c-axis) are given for illustration of possible 
magnetic structures. 
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netic structure with fc = contribute to the same neu- 
tron Bragg peaks one needs the reliable crystal structure 
data to disentangle possible magnetic contribution. For 
this purpose several data sets were collected in the single 
crystal x-ray synchrotron experiment at room tempera- 
ture and at 536K above and Tg . In addition, a limited 
slice of reciprocal space around (— |, — 1,0) satellite was 
collected at heating to identify the transition. Figure [3] 
shows a slice of the reciprocal space near [hkO] plane. The 
superstructure reflections belonging to two domains as 
shown in the figure can be easily identified. The angle be- 
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FIG. 4: Fragment of the experimental NPD pattern and the 
refined profiles for three magnetic models given by symmetry 
adapted basis functions in 14/ m space group. The full profile 
is shown by dashed blue line, the solid red lines show magnetic 
contributions to the profile. The curves for and T2 are 
shifted along y-axis for better visibility. See text for details. 



tween the domains amounted to 53.2° in accordance with 
the drawing of Fig. [2j Note that the extra peaks are cen- 
tered around forbidden nodes, because they are satelhtes 
of the Bragg peaks from the adjacent [hkl] and [/ifc — 1] 
planes. The refined structure parameters together with 
the reliability factors are given in the Table |llj Due to the 
strong absorption correction effects the atomic displace- 
ment parameters ADP can have an overall systematic 
shift. The refinement of the single crystal data in P42/n 
space group gives slightly worse, but still acceptable re- 
liability factors that we list here for completeness Rl = 
0.0860, wR2 = 0.2131 [/ > 2a{I)] and Rl = 0.1283, wR2 
= 0.2885 (all data). One can notice additional diffrac- 
tion spots at ( ^ , ^ , 0) in Fig. |3] These spots are actually 
a projection of satellite rods at C^/^J)- The in-plane 
propagation vector is equal to 5] in both average cell 
and supercell. We do not have any model to account for 
this additional superstructure, but it must correspond to 
a 2D-ordering within the (a6)-plane, e.g. an ordering of 
vacancies in Cs layers without correlations between the 
layers along c-axis. 

Using the structure data for 14/ m from the single crys- 
tal x-ray experiment we made an attempt to evaluate the 
magnetic contribution in the NPD data. We assume that 
only fully occupied Fel site has a magnetic moment. All 
the structure parameters from Table |TT] where fixed in the 
subsequent NPD refinements. Only overall ADP was in- 
troduced to account for the absorption effects in the x-ray 
experiment. The space group I4/m has eight one dimen- 
sional irreducible representations (irreps) for A; = 0, and 
all eight irreps enter three times in the magnetic repre- 




FIG. 5: Integrated intensity of (0,-2,0) superstructure satel- 
lite, lattice constants obtained from the synchrotron mea- 
surements and differential scanning calorimetry (DSC) signal 
as a function of temperature. The superstructure satellite 
intensity was obtained from the single crystal measurement 
on heating. The lattice constants were refined from powder 
diffraction synchrotron data both on heating and cooling. 



sentation for the iron in general (16i) posit ion. The irreps 
in Kovalev notatiorP^i are listed in Table lllll There are 
four complex irreps with Herring coefficient and four 
real irreps that correspond to the respective Shubnikov 
groups of 14/ m. We sorted out all the irreps and found 
that there are different magnetic configurations with the 
moment size about 2^b per iron site that do not contra- 
dict to the NPD data. The magnetic R-factors amounted 
to 17.5-24.5% for different irreps. For the illustration of 
the magnetic contribution we show in Table |III| and in 
Fig. |4] two "orthogonal" magnetic models. For T4 we 
choose the basis functions with the moments in the (ab) 
plane, whereas for r7 the moments are parallel to c-axis. 
Both models have practically the same magnetic Bragg 
R-factors 18.5 and 17.5%, respectively. We note that the 
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magnetic moment sizes on the Fel sites are not restricted 
to be the same by symmetry for complex irreps even if 
we consider the basis function only along one axis. For 
instance V't can be multiplied by an arbitrary phase fac- 
tor exp(iip) that would result in two different moment 
values. By choosing the phase (p — tt/A all the moments 
are constrained to be the same. 

The model proposed in Ref. ^ corresponds to T2 with 
the Shubnikov symbol lA/m'. Unfortunately, we do not 
observe an explicit magnetic contributions in (101) Bragg 
peak as observed in Ref. ^ in K-intercalated FeSe . This 
might be partially due to the fact that for the lattice 
constants of CsyFe2-xSe2 the (101) and (002) appear 
at the same scattering angle. In addition, in our case of 
CsyFe2_xSe2 the magnetic contribution is not so large. 
For comparison we show also the contribution of this 
model (t2, ||c) to the diffraction pattern (Fig. |4]). One 
can see that the contributions of both T7 and T2 models 
are very similar (there are small differences hardly visible 
on the figure scale), but the magnetic configurations are 
different, namely for T7 the constant moment configura- 
tion corresponding to the operators listed in the table are 
1,-1,-1,1,1,-1,-1,1, whereas for T2 1,1,1,1,-1,-1,-1,-1. Prob- 
ably the Ty model would also fit the data of Ref.l^. We 
would like to stress that the possible solutions are highly 
degenerate by the values of the R-factors and we cannot 
choose a specific model on the basis of our experimental 
data. 

Figure [5] shows the integrated intensity of the super- 
structure satellite and the lattice constants as a function 



of temperature obtained in the single crystal and powder 
synchrotron diffraction experiments. The intensity grad- 
ually disappears with transition temperature Tg ~ 500 K, 
whereas the lattice constants exhibit a pronounced hys- 
teresis indicating a first order phase transition. The unit 
cell volume is linear and does not have a visible pecu- 
liarity in the temperature region shown in Fig. [5j Inter- 
estingly, the c-lattice constant shows a decrease by 0.1% 
at the transition to the disordered phase. The crystal 
structure above is well refined in the lA/mmm model 
(Table |l]). The DSC signal has two peaks, one large at 
higher temperature and the second small one at lower 
temperature, which had been associated with the onset 
of the magnetic order from /iSR -experimenlP. The large 
DSC peak seems to be originated from the vacancy order- 
disorder transition at ~ 500 K. 
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